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Abstract 

A fundamental aspect of turbulence theory is related to the identification of realizable phase- 
space statistical descriptions able to reproduce in some suitable sense the stochastic fluid equations 
of a turbulent fluid. In particular, a major open issue is whether a purely Markovian statistical 
description of hydrodynamic turbulence actually can be achieved. Based on the formulation of a 
deterministic inverse kinetic theory (IKT) for the 3D incompressible Navier-Stokes equations, here 
we claim that such a Markovian statistical description actually exists. The approach, which in- 
volves the introduction of the local velocity probability density for the fluid (local pdf) - rather than 
the velocity-difference pdf adopted in customary approaches to homogeneous turbulence - relies 
exclusively on first principles. These include - in particular - the exact validity of the stochastic 
Navier-Stokes equations, the principle of entropy maximization and a constant H-theorem for the 
Shannon statistical entropy. As a result, the new approach affords an exact equivalence between 
Lagrangian and Eulerian formulations which permit local pdf 's which are generally non-Maxwellian 
(i.e., non-Gaussian). The theory developed is quite general and applies in principle even to turbu- 
lence regimes which are non-stationary and non-uniform in a statistical sense. 

PACS numbers: 47.10.ad,47.27.-i,05.20.Dd 
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I. INTRODUCTION 



In this and in an accompanying paper ij] the possibihty of formulating inverse kinetic 
theory (IKT) approaches to hydrodynamic turbulence theory is investigated. By definition, 
these are meant to be as phase-space models able to dehver a prescribed set (or subset) 
of fluid equations to be expressed in terms of appropriate moment equations of a suitable 
statistical equation, denoted as inverse kinetic equation (IKE). Depending of the subset of 
fluid equations to be considered different (and possibly non-unique) IKT approaches can in 
principle be developed. In the present paper, in particular, a Markovian statistical model 
of turbulence is obtained, based on the formulation of a deterministic inverse kinetic equa- 
tion for the stochastic-averaged Navier-Stokes equations. The theory permits the explicit 
construction of the local position-velocity joint probability distribution function (local pdf ) 
which advances in time the corresponding stochastic-averaged fluid fields and which can be 
shown to characterize uniquely an incompressible isothermal fluid. 

A. The unsolved problem 

The theoretical approach to the turbulence problem in incompressible fluids is one of 
the outstanding intellectual challenges of contemporary physics. Turbulence theory has 
been pioneered by Kolmogorov (K41, 2|]) who - using simple heuristic arguments based on 
dimensional analysis - first shed light on the understanding homogeneous turbulence (HT, 
i.e., stationary, spatially homogeneous and isotropic turbulence). This led to the interpre- 
tation of HT as a self-similar energy cascade in which turbulent eddies transport the kinetic 
energy of the fluid from a prescribed injection scale to a suitable dissipation scale (the so- 
called inertial range). This lead to the well-known Kolmogorov "5/3-power law" conjecture 
[2I - later confirmed by experiments and numerical results 3] - that in the inertial range 
the cascade is characterized by energy spectrum E{k) of the form E(k) = Kxo^^'^k'^/^ , 
where K^o is an universal constant called Kolmogorov's constant, k is the wavenumber, 
and n is the nonlinear cascade of energy, to be identified with the dissipation rate of the 
fluid. Subsequently the problem of homogeneous turbulence was approaches with the goal of 
quantifying the processes underlying the spectrum. Various theories were developed, for this 
purpose, to understand fluid turbulence, based on attempts to introduce suitable models for 
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the statistics of turbulent flows 0, 0]. These phenomenological theories can be developed, 
in principle, choosing for the description of fluids either Eulerian or Lagrangian viewpoints. 
The two approaches, if fluid dynamics were fully understood, should be completely equiv- 
alent. Unfortunately, at least for the treatment of turbulent fluids, we are still quite far 
from reaching this goal. The main historical reason of this situation can be understood by 
looking at the customary statistical approach based on so-called velocity probability density 
function (pdf)- method for an incompressible fluid (for a review see for example jsl). In the 
Lagrangian treatment of turbulence [denoting the so-called Lagrangian turbulence (LT)] the 
Lagrangian path R(r,t) and the velocity U(r,t) of a fluid element, initially starting at the 
position r, are determined by the equations: 

dR{r,t) 



dt 
dlJ{r,t) 



UM, (1) 
A(r,t;(7), (2) 



dt 

where the vector fleld A{r,t;g), to be suitably related to the Navier-Stokes equations, is 
assumed to depend on an appropriate statistical distribution g. Usually this is identifled with 
the joint position-velocity probability distribution (pdf) of the particle for the increments 
x(t) = R(t, y) - U(0, t)t - r, u(t) = U(t, y) - U(0, y) deflned as 

(7(x,u,t) = (5(x-x(t))(5(u-u(t))) 
/(x,u,0)=5(x)5(u). 

where the brackets denote a suitable stochastic average over a stationary statistical ensem- 
ble. For HT, the velocity increment pdf ([3]) can be shown to obey a Fokker-Planck statistical 
equation 5||. The fundamental reason why this happens is that a statistical theory of tur- 
bulence, relying solely of flrst principles, is yet not available. For this reason in the past 
various statistical models, based on heuristic assumptions about the statistics of A(t,y]g) 
and of the related pdf g, have been introduced. In turbulence theory g is usually iden- 
tifled with the velocity-difference probability density function (pdf), traditionally adopted 
for the description of homogeneous turbulence. In these cases, however, the corresponding 
Fokker-Planck equation for g generally does not deflne a closed system of moment equations 
(which should actually coincide with the Navier-Stokes equations themselves). This raises 
the critical issue of the closure of the statistical models of fluid turbulence and the difficulty 
of achieving equivalent Eulerian and Lagrangian formulations. In the past several statistical 



models have been proposed for the determination of g, which include the mapping- closure 
method, the test-function method and the field-theoretical approach . However, a critical 
aspect to turbulence theory is the possible appearance of coherent structures (like vortices 
and convective cells). In fluid turbulence the signature of the presence of coherent struc- 
tures is provided by the existence of non-Gaussian features in the probability density. During 
the last few years lots of efforts have been put into the formulation of more sophisticated 
)henomenological theories which can take into account these facts (c.f. the review article 
TJ). These approaches are based on a theoretical analysis of the Navier-Stokes equation 
or on the advanced data analysis of the exper,mentally obtamed Lagrangian path's of par- 
tides (see for example [«]). Nevertheless, despite the progress achieved in modelling key 
features of the basic phenomenology, still missing is a consistent, theory-based, statistical 
description of fluid turbulence. Clearly, such formulation - if achievable at all - should rely 
exclusively on a rigorous, deductive formulation of the turbulence-modified fluid equations 
following from the fluid equations j^. Based on a recently proposed inverse kinetic theory 
(IKT) for incompressible fluids (EUero et al, 2004-2008 [id, lul, QJ, lisj ) , here we intend to 
formulate a deterministic IKT for the stochastic-averaged Navier-Stokes equations. Key 
features of the new theory are: 1 ) the formulation is based on the local pdf, rather than the 
velocity-increment pdf usually adopted in turbulence theory (see for example [5]); 2) the pdf 
is advanced in time by means of a Markovian Vlasov-type kinetic equation; 3) the kinetic 
equation implies the exact validity of the stochastic- averaged Navier-Stokes equations (see 
below) and as a consequence the kinetic equation satisfies exact closure conditions; 4) the 
theory displays a complete equivalence of the Lagrangian and Eulerian viewpoints '[IS]; 5) 
the theory is based on first principles, i.e., besides the Navier-Stokes equations, the validity 
of the principle of entropy maximization and a constant H-theorem to be imposed on the 
Shannon statistical entropy. The theoretical setting here adopted is based on the definition 
of a "restricted" phase-space representation of the fluid [l3^, whereby the phase space (F) 
is identified with the direct-product space T = Q x V, where ^l,V ^ M^, Q is an open set 
which coincides with the fluid domain (i.e., the bounded sub-domain of M.^ occupied by the 
fluid) and V is the velocity space. 
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II. STOCHASTIC-AVERAGED FLUID EQUATIONS 



For definiteness here we consider an incompressible isothermal fluid immersed in a fluid 
domain Q (to be identifled with a bounded subdomain of with closure i7 = i7 U 5fl). 
In the set Q x I, I denoting an appropriate finite time interval / C M, the relevant fluids 
{p = Po > 0, V,p}, i.e., respectively the (constant) mass density, fluid velocity and pres- 
sure describing the fluid, are assumed to be strong solutions of the so-called incompressible 
Navier-Stokes equations (INSE) llj. The starting point of the new statistical description 
here introduced is the assumption that the fluid equations are stochastic in some sense to 
be suitably specifled. This permits us to deflne a stochastic decomposition for all the rele- 
vant quantities, i.e., the fluid flelds Z{r,t) ={V,p} as well the volume force density f(r,t). 
Namely we let, Z{r,t)= {Z{r,t)) + 6Z{r,t) and f(r,t)= (f(r,t)) + 6f{r,t), where () is a suit- 
able stochastic-averaging operator. Here, by proper deflnition of the operator () , {Z{r,t)) , 
(f (r,t)) and 6Z{r,t) and 5f(r,t) represent respectively the averaged parts and the stochastic 
fluctuations. In the sequel the precise deflnition of the operator () is not required, how- 
ever, we shall assume that it commutes with all the differential operators appearing in the 
previous fluid equations (i.e., ^, V and V^). As a consequence a suitable set of stochastic 
equations can be obtained. In particular, the equations for the stochastic-averaged flelds 
{Z{r,t)) read 

V.(V>^0. 

{N) (V) + {5N5V) = 



{stochastic-averaged INSE). Here the notation is standard [Ul]. Thus, N is the nonlinear 
Navier-Stokes differential operator A^V = + V ■ VV -|- [Vp — f] — z^V^V, where f is 
the force density, assumed to be smooth real vector fleld. and po and the kinematic viscocity 
u are real positive constants, both to be considered non-stochastic. Moreover, (A^) and 6N 
are respectively the operators (A^) (V) = g (V) ■ V (V) ^ [V {p) - (f)] - z/V^ (V) and 
SN6V = 6Y-V6Y + j-^ [V6p - 6{] - uV^6V, while Il^{r,t) = {6V6\) is the Reynolds stress 
tensor. 
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III. DETERMINISTIC IKT FORMULATION FOR THE STOCHASTIC- 
AVERAGED INSE 



The discovery of the IKT for INSE [10| suggests us to seek an inverse kinetic equation 
for the set of stochastic-averaged equations defined by Eqs.(jl]). In particular we look for a 
Markovian inverse kinetic equation (IKE) of the Vlasov-type which in Eulerian form 
reads 

L((Z))/(x,t;(Z))=0, (5) 

H{z))f = + (6) 

{Eulerian IKE). Here /(x, t; (Z)) denotes the Eulerian local pdf for Eqs.(j4]), which advances 
in time the stochastic-averaged fluid fields {Z{r,t)) and L{{Z)) is the corresponding stream- 
ing operator. We intend to construct /(x, t; (Z)) and L{{Z)) by imposing a suitable set of 
prescriptions. Besides the requirement of validity of the fluid equations (jl]), these include in 
particular: 



X = (r, v) is the state vector spanning the restricted phase space T = x V [13| and 
the vector field X has the form X((Z)) = {v, F((Z))}, where F{{Z)) = F(x, t; f, {Z)) 
is an appropriate mean-field force, to be assumed generally functionally dependent on 
the local pdf. 

By appropriate choice of the mean field force F and of the local pdf, the moment 
equations can be prescribed in such a way to satisfy identically Eqs.(jl]). For this 
purpose we assume that / is a strictly positive, suitably smooth in F x J and summable 
both in the phase-space F and in the velocity space V. In particular, we require that 
the (Shannon) entropy integral 

S{f) = - J^d^f{^,t; (Z))ln/(x,t; {Z)) (7) 

exists and that there results identically in F x / : 

(Z(r,t)) = J d'vG{^,t)f{^,t;{Z)), (8) 

where (Z(r,t)) = 1, (V(r,t)) andpi(r,t) are the velocity moments G'(x,t) = 1, v,| (u)^ . 
Moreover, u = v — V(r, t) is the relative velocity and pi(r,t) = Po{t) + {p{r,t)) is the 



kinetic pressure and -Po(^) a smooth real function (pseudo-pressure) to be suitably 
defined (see below). 

The form of the pdf / is chosen in such a way to satisfy the principle of entropy 
maximization (PEM) i.e., imposing the variational equation SS{f) = 0, with 
6'^S{f) < 0, while requiring that / belongs to a suitable functional class {/(x, t; {Z))} , 
to be determined based solely on the information available on f{'x,t;{Z)). In a 
turbulent fluid this is manifestly provided by the knowledge of the stochastic-averages 
of the fluid fields, {Z{r,t)). In this case, imposing the constraints placed by the 
moments ([8]) it is immediate to prove that PEM yields necessarily as a particular 
equilibrium solution of the inverse kinetic equation the local Maxwellian distribution 
function [kinetic equilibrium) 

fMi^,t; (Z)) = ^^exp {-X^} , (9) 

where = v^f^ = 2pi/ and (u) = v — (V(r, t)). However, in principle arbitrary 
non-Maxwellian initial kinetic distributions are possible. These are potentially rele- 
vant, in particular, for direct numerical simulations, in which the kinetic distribution 
function is simulated numerically by means of test particles. In such a case, in fact, 
small numerical errors may imply that locally the kinetic distribution function may 
actually be non-Maxwellian. 

Eq.(l5]) implies the construction of a suitable classical dynamical system, defined by a 
phase-space map 

Xo x(t) = Tt^to^o, (10) 
where is the evolution operator generated by the initial-value problem 

4x = X(x,t), 

to be viewed as the Lagrangian (or Langevin) equations for Eq. ([5]). 

The equivalence between Eulerian and Lagrangian representations. In fact the Eule- 
rian IKE can also be represented in the equivalent Lagrangian form 1^ 

J(x(t),t)/(x(t),t; {Z)) = /(x„,t„; {Z,)) = /„(x„; (Z„)), (12) 



{Lagrangian IKE) where /(x(t),t; (Z)) is the Lagrangian representation of the pdf, 
x(t) is the solution of the initial-value problem ffTTj) . /o(xo; (Zq)) is a suitably smooth 



initial pdf and J(x(t),t) 



9x(f) 
9xo 



is the the Jacobian of the map Xq x(t). 



The following theorem can be proven: 



Theorem - Markovian IKT for the stochastic-averaged INSE 

Let us assume that: Ai) Eqs.^ admit a smooth strong solution in T x I; A2) the mean-field 
force F(x,t; /; {Z)) associated to the fluid fields {Z) reads: 

F(x,t;/,(Z))=Fo + Fi, (13) 

where Fq, F^ read respectively 

Fo(x,t; /, {Z)) = - [V-n - Vpi + (f^)] + i (u . VV) + ^ (VV . u) +z/V^ (V) , (14) 

Po ^ 

F.(x,; /; (Z)) = \ („) (l^+lv • Q - J, [V-fl .q) -H |^V.fl ( J - . (15) 

where A = ^ (Po(t) + {p)) + (V) -V (Po(i) + {p)) \ A3) /(x,t; {Z)) satisfies suitable initial 
and boundary condition consistent with the initial-boundary value problem (see Ref. 

\l\]): A4) the initial pdf , /(x,to; (Zq)), is suitably smooth an strictly positive; A5) f{^,t] (Z)) 
admits the moments ([^ and ^ as well as Q = J d^vu^f and H = J d^vuuf, to be assumed 
suitably smooth; A6) the pseudo-pressure Po{t) is determined in such a way that there results 
identically in I 

= (16) 

(constant H-theorem); A7) in the fluid equations ^ and in Eqs. ( fl^ and (TR) the Reynolds 
stress tensor YIr = (5V(5V) is considered a prescribed function of {r,t) . 
It follows that: 

Bi) the velocity-moment equations of IKE ^ evaluated for the weight functions 
G{^,t) = v,|m^ coincide with the stochastic INSE equations); B2) the Maxwellian pdf ^ 
is a particular solution of IKE. In particular, the local Maxwellian distribution function 
^ is a particular solution of the IKE ^ if an only if the fluid fields { (V) , {p) } satisfy 
the stochastic-averged Eqs.^. In such a case there results identically Q = 0, n = 0; i^aj 
f{^,t; (Z)) is strictly positive in T x I. Hence it is a probability density. 
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PROOF - The proof is immediate. In fact: Bi) First, invoking Eqs. (fT3|) . (fT^ and (fTSl) . it 
follows that the velocity-moment equations of IKE Eq.([5]) for the weight functions G{r, v,t) = 



1, v,|u^ read respectively: 

V ■ (V) = 0, (17) 

I (V) + (VV ■ V) +- [V (pi) - (f)] - uV' (V) = 0, (18) 

V ■ [(V) (pi)] + (V) ■ V (pi) = (19) 



[which manifestly coincide with the fluid equations Eqs. (jlj)] . Hence, also thanks to A6 the pdf 
advances in time uniquely the stochastic-averaged fluid fields {Z{r,t)) . B2) Second,. invoking 
Al it is immediate to prove that fM{^,t; (Z)) is a particular solution of the inverse kinetic 
equation ([5]). In fact, substituting ([9]) in the inverse kinetic equation ([5]) it follows: 

L{{Z))fM = {f (V) +v • V (V)} •^/m+ 

+ {|ln(p,)+v.Vln(p,)}{g-^|/..- (20) 

-F(x, t; fM, {Z)) ■ ^-^fM + fM^ ■ F(x, t; /m, {Z)) = 0. 
Thanks to Eqs.( fT3il .(fT^ and ( |T5l) . there results: 

Li{z))fM = {I +^ ■ ^ + ;^ [V (Pi) + (f)] - 

- (u • VV) -uV' (V)} • ^^/m + /a/V ■ (V) = 

{Pl) 

which implies Eqs.dl]). Instead, if we assume that in F x /, / = /a/(x, t; {Z)) is a particular 
solution of the inverse kinetic equation, which fulfills identically the constraint equation (fT6ll . 
it follows that the fluid fields (V) , (p) are necessarily solutions of the INSE equations. If, 
instead, / = /a/ (x, t; (Z)) is a particular solution of Eq.([S]), thanks to Bi it follows that the 
stochastic INSE equations are necessarily fulfilled. B3) Finally, due to Eq.Q the entropy 
production rate reads 

|^(/W) = -^rfx^/(x,t;(Z))[l + ln/(x,t;(Z))]= (21) 
= -^rfx/(x,t;(Z))|--F(x,t;/,(Z)). 
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Hence, thanks to Eqs. (fT3ll . (fT^ . (fT5|) . the constraint Eq. (fT6i) requires 
d 



/ ^'r— — ^— - A+V • Q - - [V-nl -Q 

Jn Po{t) + {p)[ ^ p^l =i 



0, (22) 



which dehvers an ordinary differential equation for the pseudo-pressure. Assuming that the 
fluid fields (V) , (p) and the moments Q and n are suitably smooth, this equation can always 
be fulfilled in a finite time interval /. As a direct consequence, it follows that /(x,t; (Z)) is 
manifestly strictly positive in F x /. 

IV. CONCLUSIONS 

In this paper a statistical model of hydrodynamic turbulence has been formulated by 
means of a deterministic inverse kinetic theory (IKT) for the stochastic-averaged INSE. Basic 
aspects of the new theory are: A) that the IKT satisfies exactly the relevant fluid equations, 
represented by Eqs.(jl]), while B) the pdf advances in time uniquely the stochastic-averaged 
fluid fields (Z). As a main consequence, the theory fulfills identically a closure condition. 
In fact, there exists, by construction, a subset of the moment equations which is closed 
and manifestly coincides with the same set of stochastic fluid equations. The result B) 
holds provided the Reynolds stress tensor can be considered prescribed (to be identified, for 
example, with a suitably weighted velocity-space integral of the local pdf). While the form of 
the tensor remains in principle unspecified, this limitation shall be lifted in an accompanying 
paper Ij. The present theory displays several remarkable new features. In particular, unlike 
customary statistical approaches, it is based on the introduction of the local position- velocity 
joint probability density function (local pdf). In addition, the kinetic equation advancing 
in time the pdf is a Markovian Vlasov-type kinetic equation which admits a straightforward 
equivalent Lagrangian formulation. Under suitable prescriptions the form of this equation 
can be proven to be unique |12j]. Further key property is that the kinetic equation admits in 
general, besides a local Maxwellian kinetic equilibrium, also non-Gaussian solutions. In our 
view the theory provides a convenient setting both for the investigation of theoretical aspects 
of turbulence theory. These include - besides the mathematical formulation of turbulence 
problems - several important physical applications, such as, for example, the connection 
with Fokker-Planck statistical descriptions |l| and the Lagrangian dynamics of particles in 
turbulent flows and its generalization to incompressible thermofluids [l^. 
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